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Abstract 

We consider matrices containing two diagonal bands of positive entries. We show that all 
eigenvalues of such matrices are of the form r C, where r is a nonnegative real number and C is 
a pth root of unity, where p is the period of the matrix, which is computed from the distance 
between the bands. 
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Let h and k be positive integers and consider the double band matrix 



/ 





A 








• • • 

a2,fc+2 • • • 



ab+2,2 



V 



where the diagonals Ofe+j j-, are positive real numbers and the remaining entries are zero. Such 
matrices arise, for example, in second order differential equations. See [2j for an application to the 
Lame equation. The matrices as in ([T|) have a period J3, which is the greatest common divisor of 
the lengths of cycles in the directed graph r(A). Thus it is given by 



P 



5cd(6, k) ' 



'Department of Mathematics, California State University at FuUerton, FuUerton, CA 92834, USA (tm- 
cmillen@fullerton.edu) . 

^The period is also sometimes referred to as the index, or index of imprimitivity. A primitive matrix is one in 
which p — 1. 
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If b and k are relatively prime then A is irreducible. The Perron-Frobenius theorems [T] tells 
us that the spectrum of a nonnegative and irreducible matrix with period p is invariant under 
rotations by 27r/p and that there are exactly p eigenvalues with maximal modulus p{A) given by 
p{A) exp(i2irj /p), j = 0, . . . ,p — 1. The main result of this paper is to show that for matrices with 
the special form ([T|), in addition to this, the eigenvalues all lie on the lines through the pth roots 
of unity in the complex plane. 

Theorem 1.1. Let A be an n x n matrix as in ([T|). Let n = mp + q for < q < p, and let 

g = gcd{b, k). If g = 1 then A has a zero eigenvalue of multiplicity q, and mp eigenvalues 

rJ-^-ilv^ j = 0,...,p-l, s = l,...,m, (3) 

where r<j, s = 1, . . . , m, are distinct, real and positive. 

More generally, the nonzero eigenvalues of A are given by 

■ ' +1 



rs,te 



i2nj/p 



J = 0, 



,p 



h t = h 



,9, s 



1, 



n—t 



P 



(4) 



where for each fixed t, the rs,t 's are distinct, real and positive. The zero eigenvalue of A has 
multiplicity 



mod p] + {n mod g) 



[( 


n 




.9. 



+ 1 I mod p 



n 






mod p 


_9_ 



(5) 



Before proceeding to the proof we make a few remarks. First, if b = k, i.e. the bands are an 
equal distance from the diagonal, then p = 2, so all the eigenvalues are real, and come in positive 
and negative pairs. The case b = 1, k = 2 was considered in [2], where it was shown that all 
eigenvalues are of the form r exp(i27rj/3). Theorem 11.11 is a generalization of the result in that 
paper. The following corollary is an immediate consequence of the construction of the proof of the 
above theorem, and generalizes to the case when the two bands contain nonnegative and complex 
entries. 

Corollary 1.1. Let A be an n x n complex matrix with Uij = for all i — j ^ b or —k. Then the 
spectrum of A is invariant under rotations by 2tt/p. If, in addition, aij is real and nonnegative, 
then the eigenvalues of A are all of the form r e^'^'^^l'^ , where r is a nonnegative real number. 



2 Proof of Theorem 11.11 

We first prove the following lemma that allows us to reduce the problem to the case when b and k 
are relatively prime. If g > 1 we can move the diagonals "inward" at the expense of rearranging 
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them and inserting zeros. In the proofs that fohow we will employ many of the results from the 
theory of nonnegative matrices, which can be found in [3j and 



.1 



Lemma 2.1. A is cogredien^ to a direct sum of g matrices of the form ([T]) where the b and k are 
relatively prime. That is, there is a permutation matrix P such that 

f Bi • • • ^ 

B2 ■■■ 



B = PAP^ 



\ 



(6) 



where Bi is an 



1 X 



+ 1 I square matrix with non-zero entries at the positions 



(nfo + j, j) and (i, + i) and zeros everywhere else, where = h/ g and = k/g. 

Proof. We construct P explicitly. Let a be the permutation of the integers l,...,n defined as 
follows. For i = 1, . . . , define 



n — I 



a 



+ 1 



(7) 



and the partial sums A''i = 0, A^j = ni + • • • nj_i for i = 2, . . . g. A simple calculation shows that 
ni + n2 + ■ ■ ■ + Ug = n. Now, for each i = 1, . . . , g define 



aN,+j = i + {j - 1)9, j = 1,2,..., Hi 



(8) 



Then o" is a permuation of (1, ... , n), and we have 



cr.. 



k, 



(9) 
(10) 



for all pairs Ni < + j,j < Ni + Ui for some i, and likewise for Uk + i, i. Now we define the 
permutation matrix P by setting Pj,aj = 1 for j = 1, . . . ,n. Then, define B = PAP^ , so that the 
(i, j)th entry of B is bij = atr^.o-j- Since aij 7^ if and only if z — j = & or i — j = —k, as long as 
Ni < nb + j, j < Ni + Ui for some i and Ni < Uk + i,i < Ni + ni for some I, we have 



bi^n^.+i 



(11) 

(12) 



for some s and t. 



A matrix A is cogredient to B if there is a permutation matrix P such that A — PBP 



3 



Notice that each square block Ni < i,j < Ni +ni on the diagonal of B contains n/ — rn, nonzero 
elements at the positions {rib + j,j) and n/ — nonzero elements in the positions {i,nk + i), for 
a total of n — guf, = n — b elements in the positions (n^, + and n — gn^ = n — k elements in 
the positions (i, + i). Thus, since PAP^ moves each element of A to an unique position, we see 
that B has the form ([6]), as required. □ 



We can now proceed to the proof of Theorem II. 1[ 

Proof of Theorem \l.l\ Suppose, first of all, that b and k are relatively prime. Then p = b + k. 
Also, we may assume that b < k without loss of generality. Since A is irreducible with period p it 
is cogredient to a matrix in superdiagonal block form. That is, there is a permutation matrix P 
such that 



/ Ci 
C2 



C = PAP^ 



\ 





(13) 



••• Cp_i 

\Cp ••• y 

where the diagonal zero blocks are square, but the blocks Cj may be rectangular. The key to the 
proof is that we can find a P such that the matrices Cj are bidiagonal. We construct P explicitly. 
We first construct a permutation a of the integers from 1 to n. For each i = 0, . . . ,p — 1 define 

'yi,Zi, and rii by 



7i 

Hi 



ik + 1) mod p = {ik + 1) + Zi, 



V 



+ 1, 



(14) 
(15) 



and n_i = 0. An elementary calculation shows that no + • • • + rip-i = n. To simplify notation we 
also define the partial sums Ni = uq + ■ • • + Ui. Now, for each i = 0, . . . ,p — 1, define 



0-Ar,_i+j = 7i + (j - l)?*, j = 1, . . . , ni. 



(16) 



Thus a = (fJi, . . . , cj„) is a permutation of the numbers (1, . . . , n). 

Now we define the permutation matrix Pj,a{j) = 1- Then the entry of C = PAP^ is 

Cij = o-di^ay For ordered subsets a, (3 of {1, . . . ,n}, let A{a,(3) be the submatrix with rows in a 
and columns in [3. We define the submatrices Cj as follows: 



Ci 



C {{Ni.2 + 1, . . . , N,.2 + n^-i} , {iV,-i + 1, . . . , Ni_i + m}) , i = 1, . . . ,p - 1, (17) 



Cn 



C7({iVp-2 + l,...,n},{l,...,no}) 



(18) 



Note that since b and k are relatively prime, and b < k, Zi — Zj+i is either or 1. We now use this 
fact to show that each Cj is bidiagonal. 

Consider the main, lower and upper diagonals of Ci for i = 1, ... ,p — 1. Note that 



CrN,_i+j+l - CTNi+j 



li + (i - i)p - (7i+i + (i - 

-k + {zi - Zi+i)p 

-k if Zi - Zj+i = 
b ii Zi - Zj+i = 1 

li + jp - iii+i + U - '^)p) 

b + {zi- Zi+i)p 

b li Zi - Zi+i = 

b + p li Zi - Zi+i = 1 

7i + (i - i)p - iii+i + jp) 

-k + (zi - Zi+i)p-p 

-k-p ii Zi - Zj+i = 
-k if Zi - Zi+i = 1 

Now, since aij 7^ if and only Hi- j = b or i — j = —k, we see that each Ci, i = 1, . . . ,p 
upper bidiagonal if — Zj = 1 and lower bidiagonal if — Zi = 0, and that the entries 
two nonzero diagonals are all positive. 
Next, consider Cp. We have 

o"Afp-2+i - = 7p-i + (j - i)p - (1 + (i - i)p) 
= b 

<^N^-2+3 - -^J+l = 7p-l + (j - 1)P - (1 + jp) 

= -k 



(19) 



(20) 



(21) 

- 1, is 
in the 



(22) 



(23) 



Thus Cp is an upper bidiagonal np_i x no matrix. 
Now we have C in the form (jl3p . Thus 

f Di • • • ^ 
D2 ■■■ 



CP 



(24) 



\ Dp J 
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where 

Di = C1C2 • ■ ■ Cp, D2 = C2C'i ■ ■ ■ CpCi, Dp = CpCi ■ ■ ■ Cp-i (25) 

It follows that the non-zero spectra of each of the products Dj are all equal and that the non-zero 
spectra of C consists of the pth roots of the eigenvalues of Dj. Notice that each of the Dj are 
square x nj_i matrices. Furthermore, for each j, nj = m = [n/p\ or nj = m + 1, and for at 
least one j, nj = m. Thus, if m = all the eigenvalues are zero. So we assume that m > 0. 

So let j be such that rij-i = m, so that Dj is m x m. We will show that the eigenvalues of 
Dj are all real, positive and distinct. For this it is sufficient to show that Dj is oscillatory. First, 
we recall a few definitions. A matrix is totally nonnegative (positive) if all minors of all sizes are 
nonnegative (positive). A matrix X is oscillatory if there exists a positive integer s such that A"* 
is totally positive. Oscillatory matrices have the remarkable property that their eigenvalues are 
distinct, real and positive[3, cf. Ch. 6, Theorem 3.2]. Moreover, a matrix X is oscillatory if and 
only if it is (a) totally nonnegative, (b) nonsingular, and (c) satisfies Xj^j+i, Xj+i^j > for all i. We 
will verify these three conditions for Dj. 

Condition (c) follows from the fact that zq = zi = 0, so Ci is lower bidiagonal. And since Cp 
is upper bidiagonal, Dj is the product of upper and lower bidiagonal matrices. Thus, the super 
and sub-diagonal entries of Dj are all positive. Conditions (a) and (b) follow from the Cauchy- 
Binet identity. For any pair of ordered subsets a,/3 C {1,2,... ,m} of the same cardinality, the 
determinant of Dj{a,f3) is 

det {Dj{a, /?)) = det {{C.Cj^i ■ ■ ■ q[,+p-i]]) (a, /?)) 

(i+p-i) 

= E n det(q[,]](0,.i,0,)) (26) 

6»j,...,6'j+p_2 i=j 

where [[i]] = i mod p, Oj^i = a,^j+p_i = /?, and the sum is taken over all ordered subsets 
9j, . . . ,6j-^p^2 where the submatrices are defined. A simple exercise shows that all minors of 
bidiagonal matrices with positive entries on the two diagonals are nonnegative. It immediately 
follows from (j26p that Dj is totally nonnegative. Thus, condition (a) is satisfied. 

It remains to be shown that Dj is nonsingular. For this we use ()26p again, with a = P = 
{1, 2, . . . , m}. We have 

(i+p-i) 

det (Dj) = Yi det (C[[i]] (a,a)) 
i=j 

(i+p-i) 

+ E n det(q[,]](e,_i,0,)), (27) 
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where 9j-i = 0j+p_i = a. The leading principal submatrices Ci{a,a) are bidiagonal with positive 
entries on their diagonals, so the first term in (j27p is positive, and the remaining terms are 
nonnegative. Hence det (Dj) > and condition (b) is satisfied. 

Since Dj is oscillatory it has m distinct, positive eigenvalues wi, . . . , w^, and thus the spectrum 
of A consists of q zeros together with the numbers 

^i/Pe^2nj/p^ J = 1,...,P, 5 = l,...,m. (28) 

This completes the proof in the case when b and k are relatively prime. 

In the case when g = gcd(6, A;) > 1, we first form B as the direct sum of matrices as in 1^. 
The spectrum of A is the union of the spectra of the -Bj's, so we apply the previous result to each 
of the matrices Bi, . . . , Bg. Note that for i = 1, . . . ,n mod g, Bi is an ( \n/ g\ +1) x {\n/g\ +1) 
square matrix and for i = n mod g + 1, g, Bi is an [n/g\ x [n/gfj square matrix. The formulas 
dll) and dS]) are obtained by applying the result for relatively prime b and k to each of the matrices 
Bi and then counting. In the general case, although all of the nonzero eigenvalues are of the form 
^g«27rj/p^ there is no guarantee that the r's are all distinct since it may happen that the spectra of 
two or more of the Bi's overlap. □ 

In the preceeding proof the form of C = PAP^ in (I13p and ()24p depended only on the zero 
pattern of A. So, as long as Uij = outside of the two diagonal bands (6 + and (i, k + i), the 
eigenvalues of A are the pth. roots of the eigenvalues of Dj, and hence the spectrum is invariant 
under rotations by 2Tr/p. Moreover, if W6 impose the wea-ker condition that j, (X^^/jj-j-^ > 0, then 
all minors of the bidiagonal matrices Cj are still nonnegative, so the Dj^s in (j24p are all totally 
nonnegative. The eigenvalues of a totally nonnegative matrix are real and nonnegative, but not 
necessarily distinct and positive [3, cf. Ch. 6, Theorem 2.5]. Thus, we have proved Corollary ll.il 
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